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Abstract. The behaviour of a quantum rod, pivoted at its lower end on an 
impenetrable floor and restricted to moving in the vertical plane under the gravitational 
potential is studied analytically under the approximation that the rod is initially 
localized to a "small-enough" neighbourhood around the point of classical unstable 
equilibrium. It is shown that the rod evolves out of this neighbourhood. The time 
required for this to happen, i.e., the tipping time is calculated using the semi-classical 
path integral. It is shown that equilibrium is recovered in the classical limit, and that 
our calculations are consistent with the uncertainty principle. 
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1. Introduction 

It is well known that a vertical classical rod at rest, pivoted at its lower end on an 
impenetrable floor will continue to remain vertical till an external perturbation moves it 
out of its state of unstable equilibrium. Therefore, if the system were to be isolated from 
any such disturbances, the rod would never tip off. The scenario changes dramatically 
when one replaces the classical rod in the above picture with a quantum rod. In the 
latter case, the uncertainty priciple causes the rod to transit out of equilibrium, i.e., tip 
over. Many texts on quantum mechanics, for instance ask the reader to estimate 
this tipping time using the uncertainty principle; but a naive application can lead to 
incorrect results (see [2J for an interesting discussion). 

In this paper aimed at advanced undergraduates, graduate students and teachers in 
physics, we systematically analyze the evolution of a quantum rod out of unstable equi- 
librium. Our analysis involves the semi-classical approximation to path integral, and 
enables us to derive an analytical expression for the tipping time. It is also shown, 
that unstable equilibrium is recovered in the classical limit. The calculation is meant 
to explicitly illustrate and thus help students visualize, the drastic differences between 
classical and quantum mechanical evolution of systems. 

The issue of quantum mechanical rolling out of unstable equilibrium was first addressed 
in the context of an inverted pendulum in [3] . A special case of the problem we discuss 
has been considered, in the light of W.K.B approximation in [A^. A variant of the 
problem has also been addressed before, where the centre of mass of the rod is localized 
within the base of support of the rod (in our case, the base is taken to be a point), 
and the tipping time is computed numerically [21 [5]§l . However, we seek an analytical 
expression for the tipping time. 

2. Formulation of the problem 

Consider a rigid rod of uniformly distributed mass 'm' and negligible cross section, in 
contact with an impenetrable floor at the point O. We assume that the rod is constrained 
to move in the vertical plane. Let C be the centre of mass of the rod, and let the distance 
OC be 'a'. The moment of inertia / of the rod is given by / = ^ma^. The potential 
under which the particle moves is the gravitational potential given by 

V{e) = mga cos 9 9 e [-n/2, 7r/2] 
= oo 9^[-n/2,n/2] 

I In 4 , the W.K.B approximation is used to calculate energy eigenfunctions and eigenvalues, and 
the time evolved state is calculated in the limit that expectation of energy of the rod is close to the 
potential energy at the point of unstable equilibrium. 

§ In fact, the potential function used in [21 [5] is of the form Vb {cos{d — 6q) — cos(0o)), consisting of a 
minimum at ^? = between two maxima. Thus the resulting situation is a barrier tunnelling problem. 
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where 6 is measured from the vertical position. The latter condition on V^(^) follows 
from the fact that the floor is assumed to be impenetrable. The constant '(7' is the 
acceleration due to gravity. 

Since we need the rod to be localized around = initially, the wavefunction ip{6, t) at 
t = is chosen to be a Gaussian, with expectation of position < 9 >= and expectation 
of canonical momentum < / >= 0, i.e.. 



for 6 G [— 7r/2,7r/2] and zero elsewhere. Clearly, the wavefunction has the shortcoming 
that it is not continuous at 6 = ±7r/2. The justification is that the parameter a <^ 1, 
in order for the initial state to be close to classical unstable equilibrium, and therefore 
'?/'(±7r/2, 0) ~ 0. Note that with this "approximate" continuity, the wavefunction is 
normalized within the region 6 G [— 7r/2, 7r/2]. Also, we let cr be arbitrary (provided it is 
much less than one); as against [3111], where a is fixed by requiring that the expectation 
of H for ip{6, 0) is a minimum. 

Finally, we are interested in calculating the time it takes for the rod to evolve out of the 
(T-neighbourhood, to which it is localized initially. We define this as the 'tipping time'- 
ttip for the system. Note that this definition differs from the ones used before [2, HI Isfj jl. 
We identify the motion of the rod with the trajectory of the probability density maxi- 
mum in the 6-t space. This gives us a semi-classical understanding of the phenomenon 
of "tipping" of the quantum rod. The tipping time as defined here, is a measure of 
the time until which classical and quantum evolutions are close in some sense. Also, 
the definition helps us in generalising the conclusions of our calculation to all potentials 
(bounded from below) which have a point of unstable equilibrium. 

3. Semiclassical Path integral 

The time evolution operator, i.e. the propagator is defined as 



where H is the hamiltonian operator. If the eigenvalues En and eigenfunctions 0„(^) of 
H are known, then the propagator can be written as 



More importantly, the propagator is the Green function corresponding to the 
Schrodinger's time-dependent equation. Thus, knowledge of the propagator allows us 
to calculate the time evolution of any initial state under the hamiltonian H. Note also, 

II While [4] defines tipping time as the time required for the rod to fall to the floor, (21 E] use the 
average time for which the center of mass is still inside the region to which it was localized initially. 




(1) 
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that normalizability of wavefunction at all times demands unitarity of the propagator. 

There is another representation for the propagator called the path integral represen- 
tation, in which the propagator can be written as [6l [7] 

r0{t2)=e2 / • \ 

G{92,t2;9,,t,) = / D[9{t)]exp -Sm] (2) 

Je(ti)=ei V" / 

where S[9(t)] is the action corresponding to the Lagrangian L{9, 9) along the path 9{t), 



I.e. 

fi2 



S[9{t)\ = r L{9,9)dt 
Ju 



The integration in (2) is over all the paths between (6*1, ti) to (6^2 , t2)] thus the term path 
integral. It can be shown, that for ^ — >^ 0, the "classical path" which is an extremum 
for the action, provides the most dominant contribution to the path integral [6], |7j; and 
that (2) can be written as 

^(^2, t- 9,, 0) ^ G(0, t; 0, 0) exp (^^dass^ (3) 

where Sdass is the action corresponding to the classical path, and G(0, t] 0, 0) is a fac- 
tor independent of 9i and 92- In what follows, we will be using this "semi classical" 
approximation to the path integral to calculate the propagator instead of using the 
Schrodinger's equation. The justification for this will be provided towards the end of 
this section. 



Now, the fact that we need the initial state to be close to classical unstable equilib- 
rium (cr -C 1) provides for some simplifying approximations. Observe that for t < ttip, 
the wavepacket is largely contained within the a neighbourhood of the point ^ = 0. This 
suggests that for calculating ttip, one can practically replace the gravitational potential 
by y{(^) ~ mga{l — Y)for — 7i"/2 < 9 < n/2 and infinity elsewhere. 

With the above simplification, the lagrangian for the system with — 7r/2 < 9 < tt /2 
becomes 

L{9, 9) = h9^ + ^ma^u;^9^ 

where a;^ = g/a. Without the restriction on configuration space, the corresponding 
hamiltonian would clearly be unbounded from below, and thus lead to non-normalizable 
eigenfunctions. But the hard floor provides for a bounded configuration space, and thus 
one expects to maintain normalizability. 



The above lagrangian can be written in the familiar form 
L{9j) = ^Ma^9^ + lMa^n^9^ 



(4) 
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with the redefined variables M = I jo? and = (m/My/'^u. Notice that this is the 
Lagrangian for the hnear simple harmonic oscillator (SHO) analytically continued to 
imaginary frequency. This is clearly a general feature of all potentials in a small neigh- 
bourhood around a point of unstable equilibrium; a;^ being a measure of the curvature 
of the potential at that point. But the system at hand differs from SHO in that our 
configuration space is bounded. Despite this, it will be shown that one can use an SHO- 
like propagator as the time evolution operator for the given system without violating 
unitarity for "small enough" times. 

Coming back to the evaluation of the propagator, classical trajectories are expected 
to dominate the path integral as mentioned before. But because the configuration space 
is bounded in our case, there are multiple classical paths, owing to "bouncing off" 
from the boundaries. Therefore, an exact treatment would have to incorporate multi- 
instanton solutions [TJ [8]. But if t is small enough, then the direct classical path- i.e. a 
monotonic solution of the classical equation of motion 

e{r) = n^e{T) 

with the boundary conditions ^^(0) = 6i and 9{t) = 02- will dominate in comparison to 
the other classical trajectories. This path can be easily obtained by imposing the above 
boundary conditions on the general solution Ci sinh(fir) + C2 cosh(f2r), and is given by 

Oo — 9i cosh(fit) , . „ , N 
^ciass r = .1,^,, ^ sinh fir + 9i cosh fir 
smh (fit) 

The approximation clearly becomes better as t — 0. Consequently, using (3), the 
propagator can be written down as 



G{e2,t; ^1,0) ^G(0,t;0,0) 

iMa'^n 

exp 



_2hsmh{nt) 



{{el + el) cosh(fit) - 2^i^2) 



(5) 



where 



^cla.. = ((^? + ^2) COsh(fit) - 2^102) 

is the action[6j corresponding to 6'ciass given above. A complete determination of 
the propagator however, would need determination of the position independent factor 
(7(0, t; 0,0). Moreover, we need to have a clear upper bound on t for which makes 
sense. We impose unitarity on the propagator for this purpose, and claim that ([5]) is 
valid for those values of t, for which 

/ G*ie, t- e\ o)G(^', t- 0, o)d^' = b{e) 

is satisfied. But the left hand side of the above equation with the form of propagator 
suggested by ([5]), is certainly not a delta function; the reason being that not all classical 
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paths have been incorporated in Nevertheless it is plausible that one recovers a delta 
function- approximately for small enough t, and exactly in the limit t — > 0. Explicitly, 

f7r/2 

G*{e, te', o)G{e', t; o, o)de' = |G(o, t; o, o)|^ 

Tv/2 

'^''P 2;isinh(nt) ) J_^„ """P \2h^m\i{nty ) 

= |G'(0,t;0,0)|'exp {-zaO^ cosh{nt)) '^""^^^^^^^ (6) 
where a = Ma^Q/{2hsinh{Qt)). Notice that 
li,n !^ = .S(0) 

Therefore, if a — cxd, the right hand side of ([6]) approaches a delta function upto some 
multiphcative factors. Note that the phase exp (— ^6^^ cosh(fit)) goes away as 9 is set 
to zero in the exponential owing to the delta function. Finally, to recover the unitarity 
condition, one can chose G(0, t; 0, 0) appropriately (up to a phase) so as to cancel out all 
the factors multiplying the delta function. Since only |G(0,t;0,0)r will appear m our 
calculation, the phase of G{0, t; 0, 0) does not make any difference. Therefore, following 

m 

^^°'^^°'°)" U^nsmh(^^t) J 

In this sense, unitarity is respected (approximately) only for a ^ 1. This gives us an 
upper bound for t, i.e., sinh(f2t) <^ Ma'^Q/2h, for which the above results are valid. We 
state all of this in Proposition 1 below- 

Proposition 1: If t is such that sinh(fit) ^ Ma'^Q/2h, then the propagator 
G{62, t; 6i, 0) is approximately given by 



(7) 



G{e2,t; ^1,0) 



exp 



Ma'^n \ 
27cih sinh(nt) J 
iMa^n 



1/2 



((^2 + ^2^ cosh (J)t) -2^1 ^2) 



_2;isinh(fit) 

Equation ([H]) becomes exact in the limit t — 0. 



(8) 



A few comments are in order at this point, about the use of semiclassical path inte- 
gral as against the Schrodinger equation formalism. Observe that the semi-classical 
path integral allowed us to work with the approximate 9'^ potential in a small neigh- 
bourhood of the point of unstable equilibrium. This simplification is clearly not suited 
for finding the energy eigenfunctions using Schrodinger's equation. This is because in 
order to obtain the eigenfunctions, one has to impose Dirichlet boundary conditions on 
the solutions at 6* = ±7r/2, where the 9'^ approximation clearly fails. Also, one could 
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in principle solve for the eigenfunctions using the full form of the potential. The corre- 
sponding differential equation is called Mathieu 's equation (See ^ for an evaluation of 
the eigenfunctions in this manner). But the method is rather unwieldy for calculating 
the propagator. Thus, it is easier to use the semi-classical path integral, which also 
provides a more physical picture. On the other hand, our method has the limitation 
that one can only calculate small-time behaviour for systems which are well-localized 
initially. But as will be shown in the next section, this does not hamper the calculation 
of the tipping time for the rod. 

In the calculations that follow, we will use the propagator derived above in calculating 
the time evolution of the quantum rod for t much less than the required upper bound . 
This will be used to determine the tipping time, and it will be shown that our results 
are valid if the quantum rod is initially "localized enough" about the point 6 = 0. 



4. Evolution out of classical unstable equilibrium 



Proposition (1) enables us to calculate the evolution of the quantum rod whose state at 
time t = is given by the wavefunction in ([T]). The state of the system at a later time 
t is given by 

-7r/2 



n/2 



G{e,t;e',o)tij{e',o)de' 



The unitarity of the propagator ensures normalization of the wavefunction at all times 
t. For t satisfying the upper bound, one can use Proposition (1) in the above equation. 

G(0,t;0,0) 



7r(TErf(7r/2(T))V2 



7r/2 



d9' exp 



^Ma^^ ((^2 ^ ^,2) cosh(r]t) - 266') 



2a2 



(9) 



2hsmh.{^lt) 

Again, note that one is interested in the case a <^ 1. Therefore, we use the 
approximation Erf(7r/2o") ~ 1. Also, since ip{6',0) drops down very rapidly for |^| > a, 
we can extend the domain of integration to {—oo, oo). With this, the right hand side of 
([9]) becomes 
G(0,t;0,0) 



(v^a)i/2 



■ exp (ia6'^ cosh (fit)) 



d6' exp 



{ia cosh(fi)f:) 



1 



1/2 



2a2 



- 2ia66' 



(10) 



which is clearly a Gaussian integral. The integral is well defined since the coefficient of 
6'"^ term in the exponential has a negative real part, and is given by 
-G'(0,t;0,0) 



(0FfT)l/2 



exp {ia6^ cosh(fit)) 



exp 



-ia cosh(fit) + 1/2(t2)i/2 ^ [4ia cosh(fit) -2/a 



Aa'6 



2/J2 
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Having obtained ip{6,t), one can now calculate the probability density P{6,t) 



\ip{e,t)\. In doing so, we note that the phase factors in the above equation do not 
contribute. Thus, the probability density function is given by 



P(9,t) = ^\G{0,t-0,0)f 



a ' ^ ' ' ' ^' (a2 cosh' (fit) + 1/4(t4) 1/2 
exp 



1 + 4aV4 cosh^(fit). 

Remember that a also has time dependence, i.e., a = Ma'^Q/{2hsmh.{Qt)). The 
time dependence in (ITTi) can be made explicit 

Ma^na 



Pi9,t) 



^{M^a^n^a^cosh\nt) + h?smh\nt)y/^ 



exp 



(12) 

f9 smh^{nt) + M2a4fiV4 cosh2(fit) . ^ ' 

At = 0, the probability density is a monotonically decreasing function of time. For 
9^0 the function P{6,t) increases with time, from its value at t = till it reaches a 
maximum, and then decays off to zero as t oo (Figured]). The path traced by this 
maximum of probability density in the 6 — t space, can be interpreted as the classical 
motion of the centre of mass of the rod. Moreover, if 6i < 62, then the time at which 
the probability density maximum occurs at di is less than the time at which it occurs at 
6*2. Thus, the probability density maximum moves away from 6 = 0, which is what we 
interpret as the tipping of the rod. The tipping time ttip from our previous discussion, 
can be defined as the time at which the maximum occurs at = ±cr. Mathematically, 
ttip is given by 

''"^) - ° 

From 0121) and the definition above, ttip for the quantum rod is found to be 

= fi l, (n2 + M2aW)V2 ) (14) 

Note that since a ^ 1, sinh(f2ttip) ^ Ma?'VL/2h. Therefore, the hypothesis of Proposi- 
tion (1) has been shown to be true, till the rod tips. This completes the justification for 
using the propagator of Proposition (1) in the above calculation. 

We can now switch back to the original variables m, / and uj. In terms of these, 

where k = I /ma? = |. Notice that h/ma^uj ~ AdB/*^, where AdB is the reduced de 
Broglie wavelength for the rod [3]. In the quantum domain, AdB > 0. Therefore the cr^ 
term can be dropped in comparison to the other term in the denominator of (ITS!) . This 
gives 

<„„«^,i,i-(^H!!^) (16, 
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The approximately quadratic dependence on a suggests that, more locahzed is the ini- 
tial state, the faster it will tip off. It can also be checked that ttip is a monotonically 
decreasing function of u; reminiscent of the fact that stronger is the gravitational field, 
faster is the tipping. On the other hand AdB ^ cl represents the classical domain. 
Observe from (fT2l) . that in the limit ^ — > and a —>■ 0, P{6,t) for 6' 7^ 0. There- 
fore, the rod remains localized at the point 6 = 0; ergo classical equilibrium is recovered. 

It is worthwhile noticing that the above analysis would go through for any potential 
U{6) bounded from below, in the neighbourhood of a point of unstable equilibrium 60, 
with tu^ = \ U"{9o)\. Thus, provided cr ^ 1, (|T5l) and (|T6l) are generically true for all 
potentials outlined above, with slowly varying curvatures. One can also understand f|T6|) . 
at least crudely, in terms of the uncertainty principle. Note that because sinh(x) ~ x 
for X <^ 1, one can linearize (flGll as 

Since a ~ A^, the uncertainty in angular position for the initial wavefunction, and 
I [a /ttip) ~ A/, the uncertainty in angular momentum, the above equation takes the 
form 



AO.Al ^ h 
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which is in accordance with the uncertainty principle. Therefore, the tipping of the 
quantum rod can be understood as having been triggered by the uncertainty in angular 
momentum engendered by localisation of the initial state. 

5. Summary 

The quantum evolution of a rod out of unstable equilibrium was analyzed using the 
semiclassical path integral. It was shown that for small enough times, the propagator 
for the system is the same as that of a simple harmonic oscillator, analytically continued 
to imaginary frequency. This was used to compute the tipping time of the quantum 
rod. In the quantum domain, it was shown that the tipping time has an approximately 
quadratic dependence on localisation. It was also noticed that in the limit h —>■ 0, 
the classical equilibrium can be recovered, and that "tipping" in this sense is purely a 
quantum mechanical phenomenon. 
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